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ABSTRACT 

The properties of polarization in scattered light by aligned ellipsoidal grains are investigated 
with the Fredholm integral equation method (FIM) and the T- matrix method (Tmat), and the 
results are applied to the observed circular polarization in OMCl. We assume that the grains 
are composed of silicates and and ellipsoidal (oblate, prolate, or tri-axial ellipsoid) in shape with 
a typical axial ratio of 2:1. The angular dependence of circular polarization pc on directions 
of incident and scattered light is investigated with spherical harmonics and associated Legendre 
polynomials. The degree of circular polarization pc also depends on the Rayleigh reduction factor 
R which is a measure of imperfect alignment. We find that pc is approximately proportional to 
R for grains with |m|a:eg ^3 — 5, where Xeq is the dimensionless size parameter and m is the 
refractive index of the grain. Models that include those grains can explain the observed large 
circular polarization in the near-infrared, ^ 15%, in the south-east region of the BN object 
(SEBN) in OMCl, if the directions of incidence and scattering of light is optimal, and if grain 
alignment is strong, i.e., R > 0.5. Such a strong alignment cannot be explained by the Davis- 
Greenstein mechanism; we prefer instead an alternative mechanism driven by radiative torques. 
If the grains are mixed with silicates and ice, the degree of circular polarization pc decreases in 
the 3 /ini ice feature, while that of linear polarization increases. This wavelength dependence is 
different from that predicted in a process of dichroic extinction. 
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1. Introduction 

Polarimetry of young stellar objects (YSOs) 
such as T Tau and Herbig Ae/Be stars and their 
circumstellar material provides information about 
the distribution of matter and/or the configura- 
tion of magnetic field in their environment, con- 
firming that polarimetry is an important tool 
to study the physical processes in star forming 
regions. Large linear polarization up to 20% 
has been observed in many YSOs at visible and 



1 



infrar ed wavelengths (see IXamura fc Fukagawa 
2005L for a review). In linear polarization maps of 
YSOs, the polarization vectors near the star are 
often aligned in a direction parallel to the disk 
while they show a circular pattern far from the 
star. These patterns are usually interpreted as 
due to single and /or multiple scattering by cir- 
cumstellar grains jBastienlllQoi ). 

Circular polarimetry of YSOs has also been 
carried out, although fewer measurements are 
available than for linear polarimetry. Aperture 
circular polarimetry shows that circular polar- 
ization in several YSOs is small, i.e., the de- 
gree of circular polarization is » P-Ol ~ 
0. 1% in the optic a l ( Nadeau fc Bastien 1986t 



Menard, et all 119881 : iBastien et al.Ml989D . Map 



ping observations of the Chamaeleon infrared neb- 
ula (Glc dhill ct al. 1996,) and the G SS 30 reflection 
nebula ( Chrvsostomou et al. 19971 ) show circular 
polarization of 1 or 2 % in the near IR. More re- 
cently, much larger circular polarization has been 
found in near IR, i.e.,_2^ ~ 5% in the H-hand in 



R CrA (IClark et a]|l200 0h. 15% in the iC-band in 
the sout h-east region of BN (hereafter SEBN) in 
OMC l iBailev. et al.l Il998l: IChrvsostomou et aT 



2OOOI: iBuschermohle et al."2005V and 23% in the 
K-h&nd in NGC6334 (Menard,_eLal. 2000). 

At least four possible mechanisms based on 
light scattering/extinction by dust grains have 
been proposed for explaining circular polarization: 

1. Multiple scattering by non-aligned (spheri- 
cal or nonspherical) grains, 

2. dichroic extinction by aligned nonspherical 
grains, 

3. single scattering by aligned nonspherical 
grains, and 

4. multiple scattering by aligned nonspherical 
grains. 

Multiple scattering models by non-aligned grains 
(Mechanism 1) predict a circ ular polarizati on pc 
of at most a few percent ( Bastien 19961 and 
references therein), which is comparable to ob- 
served values in the Chamaeleon infrared neb- 



ula (Gledhill et all 



19961) and in the G SS 3 



reflection nebula (jChrvsostomou et al.l 119971 ). 

Therefore, one may deduce that circumstel- 
lar grains in those objects are not aligned, al- 



though low circular polarization will be expected 
even if grains are aligned, depending on vari- 
ous conditions. However, lar ger circular polar - 
ization as observed in R CrA (IClark et a ]|2000l) 



OMCl 

,2000; Busc hermohle et al 



(|Menard. et all l2000f ) 



Bai lev. et al.l 119981: IChrvsostomou et al 



20051) ■ and NGC6334 
cannot be explained by 
models without grain alignment. 

Dichroic polarization (Mechanism 2) occurs 
when light is transmitted through a medium where 
grains are aligned in a given direction. If the direc- 
tion of alignment does not change along the line 
of sight, only linear polarization is produced while 
circular polarization arises when the direction of 
alignment changes (jMartin 1978). The linear and 
circular 'interstellar polarization' observed when 
stellar light passes through d iffuse clouds is ex - 

(|2005l) 



Lucas et al 



plained by dichroic extinction, 
showed that strong circular polarization can occur 
with extinc tion, if grains arc small dielectric par- 
ticles. Buschermohle et al., (2005 ) found a correla- 
tion between J — K color and circular polarization 
Pc in the iC-band in OMCl (their Fig. 7) which led 
them to favor dichroic extinction as the mecha- 
nism for produci ng circular po l arizat ion in this 
object. However, iMinchin et al.l (jl99l[ ) suggested 
that major part of the J flux is due to scattered 
light from the Trapezium stars and free-free ra- 
diation. An inspection of 2MASS images shows 
that the J flux varies by only ?^ 0.3 mag in the re - 



gion of their Fig. 7 in iBuschermohle et al.l (2005), 
while the K flux varies by « 1 mag. Since the 
effect of extinction should be larger in J than in 
K , the J ~ K color variation does not seem to be 
due to extinction. The correlation between J — K 
and Pc can be interpreted such that the brighter 
part in the X-band, i.e., the region dominated by 
scattered light, is more polarized circularly. We 
therefore prefer the other mechanisms, i.e. Mech- 
anism 3 or 4. 

Circular polarization can also be produced by 
single scattering by aligned grains (Mechanism 3) 



single scattering by aligned grams ( iviecnanism 6 ) 
( Schmidtlll973l:lMartinlll978l:lDolginov k Mvtrophanov 



1996a; Bailev. et al.i 


1998: Chrvsostomou et al. 


2000 


Gledhill & McCalJ200r 


iMatsumura & Bastien 


2004 


Bastien & Matsumura 


20051). If the grain 



is very elongated or flattened and if the imagi- 
nary part of the refractive index is moderately 
large, circular polarization will be large, e.g. 



2 



> 30%, even for grains smaller than wavelength, 
i.e., in the R ayleigh approximation (see Fig. 5 in 



T-matrix method (Tm at) (jMishchenkol l2000t 



Mishchenko et al 



Gledhill fc McCa ll (2000)). When the size of di- 
electric grains is relatively large, circular polariza- 
tion becomes large, even for g rains which are not 
much elongated or flattened (|Gledhill fc McCal 
Those grains show different angular de- 



pendence of circular polarization fro m that for 
the Rayleigh approximation fciedhil l fc McCalll 



2000l ). Chrvsostomou et al. (|2000h explained the 



ratio of linear to circular polarization, or the el- 
lipticity of polarization, of SEEN using silicate or 
organic refractory grains with sizes of 0.1 — l.O/im. 
They ruled out metallic grains because those par- 
ticles do not explain the observed wavelength de- 
pendence of Re- 
calculations of multiple scattering by aligned 
grains (Mechanism 4) have been carried out re- 
cently (IWolf et al.l 119991: IWhitnev fc Wolfj 120021: 
Lucad l2003l) . The optical depth from the YSOs 
to the observer is usually much larger than unity, 
so multiple scattering should occur. Observations 
of thermal emission from grains in dense regions 
show significant linear polarization in the submm 
implying that those grains are nonsph erical and 
aligned ( Hildebrand fc Dragovan 1995f) . There- 
fore the study of multiple scattering (in the vis- 
ible and near IR) by aligned grains should be very 
rewarding. However these models are much more 
complex than those based on the previous three 
mechanisms. The models have now at least three 
additional parameters, two angles for the direction 
of alignment and the degree of alignment, com- 
pared to models for Mechanism 1 (multiple scat- 
tering by non-aligned grains). 

The mechanisms proposed until now to explain 
the circular polarization in SEEN are dichroic ex- 
tinction (Mechanism 2) and single scattering by 
aligned nonspherical grains (Mechanism 3). We 
know by now that multiple scattering is required 
for explainin g observations in most Y S Q env iron- 



ments fe.g.. iBastien fc MenardI Il988l . [l990h and 



therefore our ultimate goal is to study Mecha- 
nism 4. As a step in this direction, we explore 
in this paper Mechanism 3. Our results will be 
useful, among other things, for comparing results 
between single and multiple scattering, i.e.. Mech- 
anisms 3 and 4. We use two methods, the Fred- 



holm i ntegral equation me t hod (F I M) (| H olt et al 
19781 : iMatsumura fc Sekil ll99lL Il996ai bl) a^ 



2000l ). FIM can be apphed 
to tri-axial ellipsoidal particles, while Tmat is 
very efficient in evaluating the scattering prop- 
erties of spheroidal particles. We first compare 
the results of FIM with those of Tmat and show 
that the two numerical methods give essentially 
the same results under the same conditions (Sec- 
tion [271]). Since circular polarization in the pres- 
ence of weak alignment has not been investi- 
gated extensively s o far, except fo r anal ytical 
studies by [Po lginov fc Mvtrophanov ( 1978f ) and 
Dolginov fc Siklitsk^r(fr992l ). we present mod- 
els for aligned gr ains (Section 12.21). I n our 
previous papers (iMatsuniura fc Eastien 2004; 
Bastien fc Matsumura 20051 ). we showed the de- 
pendence of Pc on the scattering angle, i.e. the 
angle between the directions of incidence and scat- 
tering. Here we use spherical harmonics and as- 
sociated Legendre polynomials, and study further 
the angular dependence of circular polarization 
not only on the scattering angle but also on the 
directions of incident and scattered beams (Sec- 
tion 12. 3p . We compare our results to SEEN po- 
larization data and show that the observed linear 
and circular polarization in the K and L-bands 
can be explained if the grains are composed of 
silicates with radii of 0.15-1.5 /xm, an axial ratio 
of 2:1, and the Rayleigh reduction factor R > 0.5 
(Section l3.ip . We discuss the wavelength depen- 
dence of polarization (Section l3.2p . the grain shape 
and its degree of alignment (Section l3.3p . and the 
direction of alignment (Section 13. 4p . As a step 
toward Mechanism 4, we comment on the scat- 
tering properties if the incident beam is already 
polarized (Section l3.5p . and also on dichroic po- 
larization (Section 13. 6p . Finally, we assume that 
the grains are a mixture of silicates and ice and 
examine a possible polarization variation in the 3 
/im ice band feature (Section l3.7p . 

2. Calculations 
2.1. FIM and Tmat 

The Fredholm integral equation method (FIM) 
is one of the solutions to the light scattering 
problem, and is applicable to homogeneous and 
isotropic tri-a xial ellipsoidal particles with an axial 
ratio of a few ( Holt et al iMatsumura fc Sekil 



1991, 1996 



The light scattering process is ex- 



3 



pressed as a Fredh olm-type in t egral equation with 
a singular kernel. Holt et al. ( 1978h removed the 
singularity by using a Fourier transform which 
leads to a linear equation that can be solved nu- 
merical ly. This method i s known to be numerically 



the Stokes parameters is written as 



/. Qs Us Vs 



stable (jHolt et al.lll978r ). The major part of the 



FIM calculation is independent of the directions 
of incidence and scattering, thus FIM is efficient 
for scattering calculations in many different direc- 
tions. We have developed a numerical code using 
FIM; with our most recent version (version 2.1) 
on a desktop computer we can calculate scattering 
functions for a size parameter a;max(= 27ramax/7r) 
up to « 10, where Omax is the largest axis of the 
ellipsoid. 

A popular solution of the light scattering prob- 
lem for axisymmetrical particles is the T-matrix 



(Tma t) method ( Mishchenko 2000. : Mi shchenko et al 

Tmat expands the incident and scattered 
waves with the vector spherical wave functions, 
and their coefficients are related by a matrix called 
a " T-matrix" . One can evaluate the T-matrix 
and then solve the sca ttering p roblem numerically 
(Chap.6 in Mishch enko et all (.2000)). Although 
Tmat can be applied to particles of any shape, the 
formulation is simpler for axisymmetrical particles 
and the public domain codes are restricted to such 
particles. This is the most efficient method and 
it can be applied to particles with a size parame- 
ter up to « 100. In our simulations we used the 
Fortran program "ampld.new.f dated 04/03/2003 
written by Mishchenko. 

Both FIM and Tmat are rigorous solutions and 
give essentially the same results. As an example. 
Fig. [1] shows the degree of circular polarization pc 
(see eq.(l5]) below) for an oblate particle with an 
axial ratio of 2:1 and m — 1.7. The results from 
the two methods are in very good agreement. 

The geometrical configuration adopted in this 
paper is shown in Fig. [2l The grain is at the origin 
of the coordinate system and the direction of inci- 
dent light / is defined by the angle 9i with respect 
to the symmetry axis. The scattered light goes 
in direction S defined by {ds,(t>s)- The scattering 
angle Qsca, which is the angle between / and S", 
can be calculated by solving the spherical triangle 
AIS. 

The transformation or the Mueller matrix of 



(1) 

where j = 1, ...,4, = 1, ...,4, and the suffixes i 
and s stand for 'incident' and 'scattered', respec- 
tively. The elements Fjk are calculated with FIM 
or Tmat (see Appendix A for the sign of circular 
polarization). 

In the following calculations for oblate grains 
we first make a table of Fjk, either with FIM 
or Tmat, and then we evaluate the values of Fjk 
i n arbitrary direct ions with a Spline interpolation 
( Press et al.lll992h . The table of Fjk contains data 
of 9 X 17 X 17(= 2601) points in the 3-D param- 
eter space of (cos 6*^, cos 6*5, 0s) for 9i ^ — 90°, 
Os = 0-180°, and cfi^ = 0-180°, respectively; this 
corresponds to sampling in « 10° spacing. This 
angular resolution seems to be sufficient to obtain 
< 1% accuracy on the polarization of scattered 
light from our grain models. For prolate grains, we 
assume that the grain spins around its short axis 
and make a table of Fjk for (cos ^i, cos Os,(t>s)- For 
tri-axial ellipsoidal grains, we first evaluate Fjk 
with FIM not only for (cos 6*^, cos ^s, 0^), but also 
for the azimuthal angle (pi of the incident direc- 
tion. We then integrate Fjk for <j)i, assuming that 
the grain spins around its shortest axis and pro- 
vide a table for (cos^i, cos 6*5, (t>s)- The rest of the 
process is the same as that for oblate grains. 

We have used the following symmetry relations 
in our calculations to reduce the CPU time signifi- 
cantly. The values of F31 , F41 , F32 , E12 , ^"13 , ^23 , FiA , 
and i^24 change their sign while the others remain 
the same if the angles {9i,9s,4>s) are replaced by 
(180° - 9,, 180° - 9s,(l)s) or by (0„ 0„ 360° - 0,). 
It should be noted that the element F41 becomes 
zero if the angles {9i, 9s, (t>s) take specific values in 
the axisymmetrical models (Table 1). This implies 
that circular polarization can change sign many 
times according to the geometrical configuration. 

2.2. Imperfect Alignment 

2.2.1. Specific Cone Angle Model (S-model) 

Grain alignment would not be perfect in in- 
terstellar/circumstellar space; it may be perfect 
or imperfect, depending on alignment mechanism 
and various factors (e.g. Lazarian fc Hoanj 12007 



Hoang fc Lazarian I2OO8I for alignment by radia- 
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tive torques). To explore the effect of imperfect 
alignment, we consider two models. The first 
one is the "specific cone angle model" (S-model) 
and the second one is the "continuous distribu- 
tion of cone angle model (CD-model, discussed in 
Sect. 12.2.21 below). In the S-model, we assume 
that (1) the directions of the spin axis, the angular 
momentum, and the maximum moment of inertia 
coincide (direction A in Fig.[3|), and (2) the direc- 
tion A is distributed around another direction, the 
direction B in Fig. [3l while keeping the polar an- 
gle constant. According to assumption (1) the 
symmetry axis is parallel to A for oblate grains and 
perpendicular to it for prolate grains. For tri-axial 
ellipsoidal grains, the shortest axis coincides with 
A. Conventionally, the Rayleigh reduction factor 
R is used to express the degree of alignment and 
is related to the angle 9'^ as 



R = (Scos^e*; - l)/2 = 1 - (3/2)sin 



2 a' 



(2) 



The amount of dichroic polarization is exactly pro- 
portional to R in the Rayleigh approximation and 
is approximately proportional to R in interstellar 
grain models. Perfect alignment corresponds to 
R = 1 and non-alignment to i? = 0. 

Fig. [3] shows the relation between the various 
angles in our scattering model: the direction of 
incidence / is defined by the angle 9'^, and that of 
scattering S by the angles 9'^ and (p'^ with respect 
to the direction of alignment B. If there is no 
other scattering after this one, S corresponds to 
the direction toward the observer. The direction 
A is defined by 9'^ and (j>'^. The matrix Fjk is 
a function of 9i , 9^ , and 0s , and these angles are 
calculated with trigonometry from 9'^, 9'^, 0^, 9'^ 
and (f>'^. In the S-model, we numerically integrate 
the values Fii,...,F44 over the range (j)'^ = to 
180°. 

For natural incident Hght, i.e.. It ^ 0,Qi = 
Ui = Vi = 0, it is sufficient to consider the first 
column of the matrix Fki{k = 1,...,4). The de- 
gree of linear polarization pi, the position angle 
PA, and the degree of circular polarization pc are 
calculated with 



PA 



Pi^\jFi,+Fl/F,^ (3) 
90° - arctan(i^2i/i^3i)/2 - (4) 



where a2 is the azimuthal angle BSA in Fig l3] and 
the position angle PA is defined in the usual way, 
increasing counterclockwise as seen from the ob- 
server S with respect to the reference direction B. 
We discuss the scattering properties for polarized 
incident light in Section [3.51 

Circular polarization by pr ecessing spheroidal 
grains has been calculated by Gledhill fc McCall 
(2000), so we compare their results with ours in 
Fig. m The agreement between the two results for 
prolate grains is satisfactory since the difference 
between the two calculations is less than « 0.1%. 
However, results for oblate grains show a small 
but significant difference of « 1% whose cause is 
unknown. 

Gledhill fc McCalll (|2000l ) argued that the cir- 



cular polarization pc produced by prolate grains 
is much smaller than that by oblate grains if they 
are imperfectly aligned. They first searched for the 
direction of the maximum circular polarization in 
"perfect" alignment, i.e., the long axis is directed 
in a direction. Next they calculated pc in imper- 
fect alignment assuming the direction of maximum 
Pc is the same as that for perfect alignment. If 
the grains are relatively small, the angular depen- 
dence of Pc does not change and pc is proportional 
to the Rayleigh reduction factor R in imperfect 
alignment, as shown in Section [2.3.2l However, the 
grain size in their model (their LG model) is larger, 
and the dimensionless size parameter reaches 6 
for the maximum size in their size distribution. 
Therefore, the angular dependence of Pc varies as 
the degree of alignment changes. The prolate grain 
in the " perfect" alignment takes its maximum Pc 
value at (6*^ , 6*^ , (/.'J = (50.0°, 45.6°, 145.3°), while 
in the "perfectly spinning" alignment the maxi- 
mum occurs at {9[, 9'^,(j)'^) = (75.5°, 82.8°, 112.5°). 
The former case c orresp onds to the results of 
Gledhill fc McCalll (|200o( ) and is presented by the 



open circles and the dotted line in Fig. 21 The 
latter is shown by the broken line in Fig. |3]and is 
much larger than the former when 9'^ < 30°. Al- 
though prolate grains show lower pc than oblate 
grains, the difference is n ot as large as shown by 
Gledhill fc McCalll (|2000l ) under optimum condi- 



tion. 



Pc 



Fil/F, 



111 



(5) 
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2.2.2. Continuous Distribution of Cone Angle 
Model (CD-model) 

The S-model is simple and provides a good first 
approximation but is not sufficient to describe po- 
larization in weak alignment. In the S-model, we 
set R = 0, or 0'^ — 54.7°, for nonalignment and 
expect to find Pc — 0. However, the calculated po- 
larization is not zero but typically a few percent in 
our models. This leads us to introduce the second 
model in which the direction A is distributed not 
only in 0^, but also in 9'^. We assume that the 
cone angle 0'^ is distributed homogeneously from 
to d'^Q. Nonalignment corresponds to O'^q — 90° 
and perfect alignment appears when = 0°. It 
should be noted t hat such modeling i s not n ew and 
has been used bv IWhitnev fc Wolfj (l2002h previ- 
ously for models with — 10° and 30°. In the 
CD-model, the Rayleigh reduction factor R is re- 
lated with 6'n as 



i?=(cos2 0lo + cos0lo)/2. 



(6) 



2.3. Expansion with Spherical Harmonics 
&c Associated Legendre Polynomials 

We examine the dependence of F41 or Pci= 
F41/F11) on the angles 9i,9s, and 0s by using 
spherical harmonics and associated Legendre poly- 
nomials for perfectly aligned grains (Section l2.3.1|) 
and precessing grains (Section [2. 3. 2|) . 

2.3.1. Formulation 

If the angle 9i is constant, then the quantities 
i^4i and Pc are functions of 9a and cjjs only and can 
be written with spherical harmonics Y^{9st4>s)- 
For Fii we write 



00 n 

Y,Y.^r.raY:r{OsAs), (7) 
n— 1 m— 1 

where the coefficients hnm are calculated as 

bnm^ / F4i{9,,9s,<ps)Y;r{9s,(l>s)d9sd(l)s. 

(8) 

Since i^4i is an odd function in 0s, only odd com- 
ponents for (f)s or the terms in sin(m0s) are nec- 
essary: 



Y 



I {n - my.{2n + 1) 
27r(ri + m)l 



Pr(cos^s)sin(m0s). 



(9) 



The dependence on the angle 9i is carried by the 
coefficients bnm- We thus expand bnm with as- 
sociated Legendre polynomials P^{cos9i) where k 
is zero or a positive integer. If we choose fc = 1, 
the first term of the expansion for ^41, or the term 
with {n,m,l) — (1, 1, 2), coincides with the expres- 
sion for the Rayleigh approximation. We therefore 
set k = 1, and write bnm as 



b: 

where 



1=1 



21 + 1 

^anmiJ ^^^^^^^ Pl [cos 9,), (10) 



(11) 



21 + 1 

annii = J ^j^^j^ _^ -^^ J bnm{9,) Pi {cos 9i)d9, 
Combining eqs.(7)-(ll) yields 

00 00 n 

F4i(a„ 0s, 0s) = ^ ^ ^ anrnir'^'ids, c|)s,e^), 

(12) 



l—l n—1 m—1 



where 

anml — / / / Fil{(^i^Qs,4>s) 

J 01=0 J0s=O J 9^=0 

a 



y.r""\0s,<Ps,e,)d9sdcj>sd9,, (13) 



and 

rnml ( 



{21 + l){n - m.)\{2n + 1) 
Y AtiI{1 + l){n + m)\ 

X (cos (cos 61s )sin(m0s). (14) 

The functions J"™' as given by ea. p^ are or- 
thonormal, 



7r />27r />7r 



/""'(&„ 0s, 0s) 



.0 J4>s^o JOs^o 

1 m'l' 



x/" ™ ' (6*,, ^s, (t)s)d9sd(t)sd9, = 5„„'(5w'5h', 

(15) 

where 5nn' is Kronecker's ^-function. 
With this relation 

P;"(-cos0s)P,'(-cos0,) = 

-(-l)"+™+'Pr(cos0s)Pi'(cos0,), (16) 
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we see that the product {cos 9 s)Pl {cos 9 i) in 
eq. (|14p changes sign when n + m + I is even but 
not when it is odd. The quantities pc and Fn show 
the same properties to the case oi n + m + I even, 
i.e., they change sign when 0i is replaced hy tt — 9i 
and 9s hy n — 9s, while (j)s unchanged. Therefore 
the coefficients anmi vanish when n + m + l is odd, 
and Pc and F41 can be written with only the even 
terms in n + m + I. As a result we can write the 
first few terms of F41 as 

Fii{9,,9s,^s) = 

3 [5 

- -i / — aii2 sin 9i cos 9i sin 9s sin (hs + 

4 V TT 

3 [5 



{17) 



— \ —0211 sm Ui sm Us cos f ^ sm <ps + 
4 V TT 

f 5 [Y 

^\ —0.222 sin 0i cos 9i sin^ 6*^ sin 2(?!)s + 

8 V TT 



When the Rayleigh approximation is valid, only 
the first term (0112) is necessary for the expression 

of i^41. 

The degree of circular polarization pc can be 
expressed the same way as explained above and 
many examples are shown in Fig. [5] For dielectric 
particles, the overall dependence of Unmi on \m\xeq 
is similar, except for two cases: prolate grains com- 
puted with the "perfect" alignment (Fig. [Sji) and 
oblate grains with axial ratio of 3:1 (Fig.[S^). The 
term in aii2 is the dominant one for |m|a:eg ^ 2, 
which indicates that the Rayleigh approximation 
holds. Around jmlxg^ ~ 2.6, the sign of 0112 
changes, and for jmlxeg > 2.6, both the terms in 
0211 and aii2 are large. Since \m\xeq is a size pa- 
rameter normalized by the wavelength inside the 
particles, the properties of circular polarization 
are determined by the transmitted wave within 
the particles. Absolute values of 0112 and 0211 
are largest when |m|a;eg « 3 or 4 and they start to 
decrease for larger values of Xeq, for the cases we 
computed for dielectric particles with m = 1.7 or 
1.7 + O.li. 

The values of a„,„; for metallic particles are 
functions of Xeq, and not of mXeq, i.e. the overall 
dependence of their a„,„/ on Xeq is similar (Fig.[5).i 
and i). Therefore, the mechanism to produce cir- 
cular polarization in metallic particles seems dif- 
ferent from that in the dielectric ones. The abso- 
lute values of 0112 and 0211 are largest for x^q ~ 1 



and they decrease when Xeq ^ 1 or x^q ~ 0. 

2.3.2. Precessing Grains 

When the grain is incompletely aligned, the val- 
ues of F41 and Pc decrease, and the angular depen- 
dencies change. For precessing grains (S-model, 
Section 12.2. ip , we change the Rayleigh reduction 
factor R, and investigate the variation of the func- 
tion /"'"'(6'^,0'^,6'O (eq.dlll)). We write the mod- 
ified functions as g"-"^'^{9'^,9's,(j)'s,R) under the ef- 
fect of precession, which is identical to J"™' for 
R = 1. We expand (^"™' with spherical harmonics 
and associated Legendre polynomials already in- 
cluded in J"™' J in the same way as explained in 
Section 



00 00 



I' — l n' — l m' — l 



(18) 



For {n,m,l) = (1,1,2) or (2,1,1), we find that 
gumi jjg^j^ written as 

g"™'(0^, 9's, i?) - R ■ r"'\9[X, ^'s)- (19) 



The angular dependence in these cases does not 
change and the amplitude is proportional to R un- 
der incomplete alignment. 

Terms with higher (n, m, I) are required for the 
expansion of g^™^ for higher values of (n, m, I), i.e. 
{n,m,l) ^ (1,1,2) nor (2,1,1). The examples of 
aJ^J^,;, (i?) shown in Fig. [6] have been numerically 
evaluated for {n,m,l) = (2,2,2), (2,1,3), (3,1,2), 
(3,2,1), (3,2,3), and (3,3,2). As expected, the 
value of a"™; is unity for perfect alignment {R = 
1). The values of a^^j;^,;, for {n' , m' , I') ^ (n, m, 
are not large especially when R is small. This 
means that large circular polarization cannot be 
expected from large grains if alignment is poor. 

2.4. Results for Imperfectly Aligned Grains 

To evaluate how the results change when 
grains are not perfectly aligned, we assume the 
S-modcl (Section 12.2. ip and calculate the means 
of F4i{R)/F4i{l) (Figs. [7^ and c) and the stan- 
dard deviations (Figs. [7|3 and d) for values of 
R = 0.25,0.5 or 0.75 for oblate grains and one 
case for prolate grains. We assume that the sym- 
metry axis of prolate grains is perpendicular to 
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spin axis, and calculate the ratios of values F41 
for those grains over F41 in perfect alignment. 
This condition is expressed as i? = —0.5 in Fig. [71 
We took 1785 cases or less into account out of 
2601 (see Section 12. ip because we excluded 816 
cases where -F4i(l) — due to symmetry (see Ta- 
ble 1) and also other cases where ^4i(l) ~ 0. In 
the Rayleigh approximation or when mxeq ^ 2, 
the means of i^4i(i?)/i^4i(l) are nearly equal to R 
and the standard deviations are relatively small, 
as expected (Section [2.3.2p . For larger mXeq, the 
means become smaller, and the standard devia- 
tions larger. However, the means for m — 1.7 



are nearly equal to R up to mxe 



4, except 



for a drop near mXeq ~ 2.6 (Fig. [7^). The value 
of F41 for m — 1.7 can be expressed well by us- 
ing only the component of 0112 and 0211 (FiglS^). 
Since these components are proportional to R, 
F41 is also proportional to R. At mXeq ~ 2.6, this 
proportionality fails because aiu is zero and the 
relative contributions from other terms become 
large. The maximum circular polarization \pc\ is 
large when mXeq ^ 3 while the degree of linear 
polarization Pi is small (FigslTfc and f). 

For the dielectric particles considered here, one 
or two peak(s) is (are) expected for pc as a function 
of the angles {9i, 9s,4>s), because the values of aii2 
and 0211 are much larger than higher order terms, 
as shown above. It is highly possible that one of 
the peaks in the Pc can explain the high circular 
polarization observed in SEEN of OMCl. We thus 
investigate the directions which correspond to the 
maximum value of pc as a function of mXeq for 
particles with an axial ratio of 2:1, in a parameter 
space of {0^,0,, = (0-180°, 0-180°, 0-180°) 
(Fig. [5]). The values of angles for minimum pc, i-e. 
the negative of maximum value, can be retrieved 
by using a symmetry relation (Section l2.ip . In the 
Rayleigh approximation, the circular polarization 
Pc is expected to reach its maximum value when 
i0^,0s,(t)s) = (135°, 90°, 90°) (see eq.(Bl) in Ap- 
pendix B), while Figs. [Hb-d show (129°, 90°, 90°) 
when mXeq ~ 0. The difference in the values 
of 0i is due to the coarse sampling in the angle. 
One can see systematic variations of the angles 
corresponding to maximum pc value as a func- 
tion of mxeq, although those variations are not 
clear for larger mXeq. For mXeq « 3 — 5, the 
value of Pc takes its maximum in {0i, 0^,(1)3) ~ 
(40 - 60°, 40 - 90°, 100 - 140°), and this property 



does not depend on shape considered here. 

3. Discussion: the SEEN Region in OMCl 

We now compare our results to the observed 
large polarization in SEEN o f OMCl. We use the 



circu lar polarimetry data by IChrvsostomou et al.l 
(I20p0 | l (their Tabl e 2), a nd the linear polarimetry 
bv iMinchin et al.l (|l99ll ) (their Table 4). These 
tabulated values are observed in the box region of 
size 3" X 3" centered at 22" east and 6" south from 
the EN object. 

3.1. Grain Parameters 

A few papers have studied the parameters of 
grains which can explain the observed polarization 



i n the SEEN region so far. IChrvsostomou et al 
(l2000[ ) showed that oblate grains of silicate with 
an axial ratio of 2:1 can explain the ratio of circular 
to linear polarization observed in SEEN provided 
the lower cut-off and the upper cut-off in a power 
law size distribution are set to amin — 0.1/im 
and Qmax = 1.0/im, where the size is the equiv- 
alent radius of a sphere with the same volume. 
Lucas et al. I (I2OO5I) found that a distribution of sil- 
icate grains with sizes up to amax — 0.75/im could 
produce circular polarization in the K-hand if 
their axial ratio is larg er than 3:1. However , repro 



ducing the model by Chrvsostomou et al. ( 2000h 



or that by iLucas et al.l ( 20051 ). we find that the 



degree of circular polarization pc is too small to 
explain the observed one {pc = 9.4%) in the L- 
band (A = 3.6/xm), though the circular polariza- 
tion Pc = 15% in the X-band can be explained. 
The angular dependence of those grains varies so 
greatly from the K to the L-band that it is impos- 
sible to explain the observed circular polarization 
in the K and L-bands simultaneously if the an- 
gles of incidence and scattering do not depend on 
wavelength. 

We thus examine the effects of grain size pa- 
rameters in a power law size distribution and have 
fomid that the upper cut-off amax is the most effec- 
tive for changing the circular polarization pc in the 
L-band, while amin and the power index are less 
effective. Fig. [S] shows the circular polarization pc 
of oblate grains with axial ratios of 2:1 or 3:1 as a 
function of the upper cut-off a^nax , setting amin as 
o-min — O.lamax, and the power index as —3.5. If 
the axial ratio is 2:1 (solid lines), the observed val- 
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ues of polarization {pi = 38% and Pc — 15%) in the 
JsT-band can be explained if a„iax > 1.0/im, while 
those in the L-band {pi — 57% andpc = 9.4%) can 
be explained only if Umax ^ 1.5/im. The grains 
with the axial ratio 3:1 (dashed lines) can explain 
the observed circular polarization in the L-band 
even if amax ~ 1/im. However, the calculated lin- 
ear polarization is higher than the observations 
in both the K and L-bands. Therefore, we set 
dmax = 1-5/xm, amin = O.lS/Ltm, and the axial ra- 
tio is 2:1 as a typical grain model in the following 
calculations. 

As for the refractive index, we use the data for 

" smoothed astronomical silicate" ('Wein gartner fc Draind 
[2OOI: Draine.2003,) . Although the refractive index 
m of the silicate depends little on wavelength A 
in 0.5 < A < 5/im, eg. m = 1.661 -t- 0.035i in the 
JsT-band (A = 2.2/^m) and m = 1.638 + OMli in 
the L-band (A = 3.6/im), the "effective" refrac- 
tive index shows large variation around A ~ 3/im 
if ice is mixed in. We investigate this effect in 
Section 13.71 



3.2. Wavelength Dependence 

We show the wavelength dependence oipi, PA, 
Pc, and Pc/pi in Fig. [TOl for models of oblate, pro- 
late, and ellipsoidal grains composed of silicate 
with an axial ratio of 2:1 {2:^/2:1 for ellipsoids). 
Here we use the CD-model fSection l2.2.2p for grain 
alignment, and we take the Rayleigh reduction fac- 
tor R as 0.5 or 1.0. We assume optimum directions 
of incidence and scattering to obtain the largest 
circular polarization. The overall trends of the ob- 
served quantities seem to be explained. However, 
the observed linear polarization pi is systemati- 
cally larger than the computed values in most of 
our models. If we decrease the size, the fit for pi 
becomes better, but the circular polarization pc 
is smaller than the observed values. The agree- 
ment for position angle PA is not good (Fig. [TUb) 
and the reason is not clear. The values of both pi 
and Pc drop in the iJ-band, compared with those 
in the K-hand. Our models cannot explain such 
sudden variations, as the scattering properties of 
dielectric particles are less sensitive to wavelength. 
Since metallic particles show more dependence on 
the size parameter Xcq (Figs.[5h and i), such par- 
ticles may explain the observation if they exist 
and are aligned. An alternative possibility would 
be contamination by natural light in the i?-band. 



Minchin et al.l (|l99ll ) pointed out that scattered 
light from Trapezium stars and free-free radiation 
dominate over the scattered light from IRc2 in the 
J-band. It is possible that a similar situation may 
occur also in the H-hand to some extent. 

3.3. Degree of Alignment 

We explore the range of the Rayleigh reduction 
factor R which can explain the linear and circular 
polarization observed in the K and L-bands for 
the CD-model. Using the results of calculations 
in Section H we examine 4913 (= 17 x 17 x 17) 
models in the 3-D parameter space of 6i = 0—180°, 
0^ = 0-180°, and (/f^ = 0-180°. Assuming a value 
for the Rayleigh reduction factor R (column 1 in 
Table [5] and [3]) , we count the numbers out of 4913 
models which satisfy the following conditions: 



Pc 



Pi 



5pc 



, ^ obs 
<Pc<Pc 



Spc, 



and 



dpi <pi Kpf" +6pi. 



^p^obs _ < SPA, 



(20) 
(21) 



(22) 



and the results are shown in column 3 in those 
tables. In the X-band (Tabled, the observed val- 
ues (and adopted ranges of acceptance within the 
parentheses) are " = 15(3)%, pf" = 38(7)%, 
and PA°'" = 0(5)°, and in the L-band (Tabled, 
those values are pf" = 9.4(3)%, pf' = 57(11)%, 
and PA°'"' = 0(5)°. Those adopted ranges of ac- 
ceptance are larger than the observational errors, 
because the latter are too small to yield meaning- 
ful results with our models. The column 2 in Ta- 
ble [5] and [3] is the maximum \pc\, and the columns 
4, 5, and 6 are the averages (and standard devia- 
tions within the parentheses) of the angles 9i, 9s, 
and (j)s, respectively, with which the model satis- 
fies the above conditions. 

Tables [D and [3] show that R should be larger 
than ^0.4 or 0.6 to explain the observations. 
These values are much larger than the value of i? w 
0.25 which is de rived for oblate grains w i th th e 
axial ratio 1.5 bv iHildebrand &: DragovanI (|l995r ). 
from dichroic linear polarization at A = 2.2/xm, 
thermal emission at 100/im, and also from the ob- 
served 9.7 fim silicate absorption band feature. It 
is also noted that the fraction of models which sat- 
isfy the observations is small, i.e., « 10 to 80 out 
of 4913, or w 10-3 to 10"^^ This fraction will 
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increase if we adopt larger values for the ranges 
of acceptance, though it is still « 10~^ if the 
ranges are two times larger than the present as- 
sumption. This means that the optimal directions 
of incidence and scattering of light are required to 
explain the observations, and that large circular 
polarization would not always be expected even if 
the alignment is strong. Alternatively, it would 
be necessary to include other light scattering pro- 
cesses that are not considered here, i.e., multiple 
scattering and/or dichroic extinction, to fully ex- 
plain the observations. The effects of those pro- 
cesses are discussed in Sections 13.51 and 13.61 

Fig. [Tib shows the Rayleigh reduction factor R 
necessary to obtain pc = 15% in the X-band as a 
function of the axial ratio r, where the axial ratios 
for ellipsoids are set as r : y/r : 1. We show exact 
results computed with the CD-model (lines and 
open symbols). Also shown are results from the 
approximate expression R = 0.15/^™"^ (crosses 
and pluses), which is based on the fact that circu- 
lar polarization is approximately proportional to 
R if mxeq < 5 (Section [2X2]) . The Rayleigh re- 
duction factor R shows a rapid decrease for r « 1, 
reaches a minimum around r w 2, and then slowly 
increases with r for r > 2. This figure shows again 
that strong alignment, i.e., R > 0.5, is required to 
explain the observation, even if the shape is highly 
elongated or flattened. 

How can we explain such large values of R? We 
give a few comments here. The Davis-Greenstein 
mechanism, which is based on grain rotation 
driven by gas-grain collisions and on paramag- 
netic relaxation, gives a moderately large value of 
the Rayleigh reduction factor R when the ratio of 
grain temperature to gas te mperature is signifi- 
cantl y small, i.e., << 10^^, (jRoberge fc Lazarian 
19991 ). However, in dense environments such as 
the OMCl region, the gas and grain temperatures 
should be nearly equal. Thus we cannot expect R 
to be large with the Davis-Greenstein mechanism. 
The alignment by the difference of velocities be- 
tween gas and grain, originally known as the Gold- 
mechanism, is expected in the presence of A lfvenic 
or magnetosonic waves (jLazarian I1994L [l997.1 . 
The maximum R expected with this mechanism 
is « 0.25 and is smaller than what our analysis 
yields. The alignment by radiative torques is ex- 
pected when irr egularly, or helically, shaped grains 
scatter light (jDolginov fc: Mvtrophanov 



Draine &: Weingartner 119961: Lazarian fc Hoaii3 
2007 ; Hoang fc Lazariaiil 2008|). For a revie w of re- 
cent alignment theory, see LazarianI ( 2007 ). Since 
SEEN is a reflection nebula, the grains in this 
region are strongly illuminated by the star IRc2. 
Thus we may expect this mechanism to be work- 
ing there. It is unfortunate that the study of 
this mechanism is beyond the scope of the present 
paper because we cannot evaluate the effects pro- 
duced by helical grains with the FIM nor the Tmat 
method. 

3.4. Direction of the Alignment 

We have discussed the angles of incidence {9[) 
and scattering {9'^ and with respect to the ref- 
erence frame of the grain. However, in the astro- 
nomical context, it is much more convenient to use 
other angles, i.e., the position angle PA' that is 
an azimuthal angle of the alignment projected on 
sky, measured from north to east, and the scatter- 
ing angle Qsca- We need one more angle to speciiy 
the direction of alignment in space, i.e. the incli- 
nation angle between the alignment and the line 
of sight, and this angle is the same as 9s already 
defined above. We assume that the SEEN region 
is illuminated by IRc2 from a direction with po- 
sition angle 84°, i.e. illuminated almost from the 
west direction. 

With trigonometry, we calculate the angles PA' 
and Qsca in columns 7 and 8, respectively, in Ta- 
bles [2] (the if-band) and [3] (the i-band). Almost 
all the models with an axial ratio of 2:1 and the 
oblate model with 3:1 show similar results if their 
size distribution is aeq = 0.15 — 1.5/Ltm. How- 
ever, the results of prolate with 3:1 are differ- 
ent from the others, suggesting a strong shape 
dependence for elongated particles. The models 
with smaller size distribution of aeq = 0.1 — 1/im 
also show different results, especially in the L- 
band. With the exception of these models, the an- 
gles of acceptable models are {9'^, 0s) ~ (60° — 
80°, 40° - 60°, 120° - 140°), PA' « 20° - 50°, 
and esca « 100° - 110°. The deduced posi- 
tion angle PA' « 20° — 50° is almost perpen- 
dicular to the direction of the magnetic field, i.e. 



PA' 



150°, observed in the region of SEEN 



([Chrvsostomou et al.Ml994l ). and also to the av 
erage dire ction, PA' ft! 120°, in the overall region 
of OM Cl (jChrvsostomou et al.iri994l : lHoude et al ' 



20041 ) ■ This suggests that the magnetic field would 
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vary locally in a small region of SEBN where light 
is scattered, if the direction of alignment is parallel 
to the magnetic field as is usually assumed. 

3.5. Effects of Polarized Incident Light 

We now assume that the incident light is 
polarized, and investigate how the polarization 
status is changed. The normalized Stokes pa- 
rameters {q'g,u'g,v's), where q'^ = Q's/I's etc., 
of scattered light for polarized incident light 
are calculated with eq.ll]), and are compared 
with those {qs,Us,Vs) for nonpolarized incident 
light. Figs. [T^a) and (b) show the results for 
{qi,Ui,v,) = (0.3,0,0), Figs, injc) and (d) for 
(0,0.3,0), and Figs, [lite) and (f) for (0,0,0.3). 
We assume oblate grains with an axial ratio of 
2:1, composed of silicate, and size distribution of 
aeq = 0.15 — 1.5^m. While the Stokes parameters 
q'g and u'^ show significant deviations, w30%, from 
qs and Ug (Figs. fl2k and c), the effect for circular 
polarization v'^ is less significant, «10% (Figs.fT2b 
and d). Since the intensity of scattered light 
is approximated by « -Fii/i, the expressions of 
the scattered light become simpler and are given 

by 

q's ~ /2I + .f22qi + f23Ui + f24Vi, (23) 
u's ~ /3I + f32qi + f33Ui + f34Vi, (24) 



and 



/41 + f42qi + fi3Ui + fiiVi, (25) 



where fjk — Fjk/Fn. Since fjk varies from -1 
to 1, we can set upper and lower limits for the 
deviations, and show them as the dashed lines in 
Fig. [121 The quantities q'^ and u'^ go to the limits, 
while the values of v'^ are far smaller. Therefore, if 
the incident circular polarization Vi ~ 0, we may 
write the deviation of circular polarization from 
the nonpolarized incident model as 



(26) 



where pu {= ^/qf+uf) is the degree of linear po- 
larization of the incident light. The effect of a 
polarized incident light is not significant unless pu 
is extremely large. 

Fig. [131 shows contour plots of the number of 
models that satisfy the observed linear and/or cir- 
cular polarization in the if-band (eqs.(20)-(22) in 
Section [33]) out of 4913 models (see Section IXSj) . 



for oblate grains with an axial ratio 2:1, com- 
posed of silicate, with size distribution of = 
0.15 — 1.5/^m. Here we assume that the normal- 
ized Stokes parameters of incident light qi and Ui 
go from —1 to 1, and that the circular polariza- 
tion Vi is null. The number of models with larger 
linear (Fig. [T3k ) and circular (Fig. [T3h ) polariza- 
tion increases with larger qi and Ui, while the 
models that explain the position angle PA are re- 
stricted only to a region of relatively small qi and 
Ui (Fig. [T3b). Therefore, the acceptable models 
that explain all the observed properties are also 
found in a region of small qi and u^, centered on 
(q^Ui) = (0.15,-0.27) (Fig.[T3li). In this region, 
the number of acceptable models is about 50, and 
is two times larger than that on {qi,Ui) = (0,0), 
showing that the presence of optimum polarized 
incidence is effective to explain the observations. 
However, it is also noted that very large incident 
linear polarization, i.e. > 50%, is not useful. Al- 
though the effects of polarized incident light are 
significant, large incident linear polarization will 
strongly affect the position angles, not consistent 
with the observation. 

3.6. Polarization by Dichroic Extinction 

If grains are strongly aligned as discussed in 
Section 13.31 the transmitted light through space 
containing those grains would be also linearly po- 
larized. We estimate the ratio of dichroic linear 
polarization to optical depth p/r in our models, 
as a function of the axial ratio r of the grains 
(Fig. [Tib) . We show results obtained with formula 
p/t — R - {p/T)per, where R is the Rayleigh reduc- 
tion factor necessary to obtain pc = 15% in the 
K-hand as given by our exact calculations, and 
{p/T)per is the value for perfect alignment. The 
value of {p/T)per becomes larger with r, and thus 
the derived value oip/r also becomes larger, even 
if R decreases. When r w 2 — 2.5, the value oip/r 
is nearly equal t o the observed maximum w 0.07 
(Fig.l of Jones! ( 1989I )). although the exact value 
of p/t in SEBN is not known. 

Circular polarization can be produced by ex- 
tinction with aligned nonspherical grains, i.e., 
dichroic polarization or Mechanism 2 explained 
in Section 1. The exact solution of the Stokes 
parameters with this mech anism is found i n 
Whitnev fc W)ii (|2002h and iLiicas etldl (l2005l ). 
We can write an expression of circular polarization 
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as 



pi^u'TRKs^/Kn, 



(27) 



in the first-order approximation, where is a 
normahzed Stokes parameter " m" of the incidence 
light, T is the optical depth along the line of 
sight, R is the Rayleigh reduction factor, and Ku 
and Km are components of the extinction matrix 
(jMishchenkG et al.ll200d ). The value of K34/K11 
is « 0.1 in our spheroidal/ellipsoidal sihcate grains 
with radii of 0.15-1.5 /xm at A = 2.2/im. If we fur- 
ther assume « 0.3, r w 1, and R « 0.5, then 
we obtain pj? ft! 1.5%. Although each parameter is 
quite uncertain, the calculated value in our model 
is far below the observed circular polarization of 
15%. 

If the size of grains is much smaller than wave- 
length, i.e. within the Rayleigh approximation, 
the effect of dichroic polarization will increase, in 
contrast to scattering. An example is the model 
of those grains that are assumed in diffuse inter- 
stellar space, i.e., Qmax ~ 0.25/im. As for lin- 
ear polarization, the value of (j>/T)per will be 0.28 
for oblate silicate grains with an axial ratio of 
2, and the observed maximur n valu e of p/r « 
0.07 in various clouds ( J one^ Il989l ) can be ex- 
plained with R K, 0.25. For circular polarization, 
the va lue of KM|K^^ w ill be much larger than 



unity (jLucas et al 



200a|). For dielectric particles, 

4 



is proportional to Oeg, while Knvariesasa^^, 

goes as in the 



and thus the ratio K^i^^/Kn 
Rayleigh approximation. If dielectric small non- 
spherical grains are aligned, then the value of 
1^34 /ifii would be larger than unity, and the ob- 
servations could be explained with dichroic polar- 
ization. Therefore, on the assumption of smaller 
grains, one can obtain large linear and circular po- 
lariz ation with the mech anism of dichroic polariza- 
tion (jLucas et al.ll2005l ). 



3.7. The 3/im Ice-band Feature 

Spectropolarimetry of transmitted light of BN 
object and other IR sourc es shows that the S/xm ice 
band is linearly po larized ( Dick fc Beichma"nlll974l : 
iHough et al. 1996f l. This polarization is explained 
with the accretion of ice on aligned spheroidal si l- 
icate and / or graphite grains (|Lee fc Drainelll985h . 
We assume that such grains with ice are present 
also in the SEEN region, and examine the effects 
of the 3/xm ice band on polarization in scattered 



light. In our model, the grains grow in size by 
accretion of ice, and the refractive index is as- 
sumed homogeneous and is approximated with the 
Bruggesman mixing rule which is bas ed on an ef- 
fective medium theory (jKrtigel l2003h . This ap- 
proximation is adopted because our present codes 
of FIM or Tmat cannot calculate ellipsoidal core- 
mantle grains. We use the refractive index for 
the " strong ice mixture" w ith temperature 120K 
from lHudgins. et al. I (|l993h . which contains H2O : 
CH3OH : CO : NH3 = 100 : 50 : 1 : 1. The 
refractive index for silicate is the same as in the 
previous sections. 

Fig. [14] shows the results for oblate grains with 
an axial ratio 2:1. The size of silicate grains with- 
out ice is 0.15 — 1.5/im. The degree of linear po- 
larization pi in the 3/im ice feature increases with 
the volume of ice, while the circular polarization 
Pr. decreases. In t he dichroic extinction model 
( Aitken. et al. 20061 ). both the degrees ofpi and Pc 
increase in the feature. The difference is due not 
only to the light scattering process itself, but also 
to the difference of grain size in the models, i.e. 
grains are larger in our model. Therefore circular 
polarization in the 3/im band is a useful diagnos- 
tic for finding out which process is working in the 
SEBN region. 

4. Conclusions 

We have studied polarization in scattered light 
by imperfectly aligned spheroidal or ellipsoidal 
grains with FIM and Tmat. The shapes of grains 
considered here include tri-axial ellipsoid which 
has not been investigated before. Our main con- 
clusions are as follows: 

1. With using spherical harmonics and associ- 
ated Legendre polynomials, we have investigated 
the angular dependence of circular polarization pc 
or the component F41 of Mueller scattering ma- 
trix. For dielectric grains that are aligned spinning 
around shortest axis, the angular dependence of Pc 
or F41 does not vary much in different shapes, if 
the grains are not much elongated and not large, 
i.e., the axial ratio of grains is up to about 2:1 
and \m\xeq is up to » 5. For those grains, Pc 
is approximately proportional to the Rayleigh re- 
duction factor R, even when the scattering proper- 
ties are far from the Rayleigh approximation, i.e. 

\m\Xeq « 3 - 5. 
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2. To explain the large linear and circular po- The refractive index of smoothed astronomical 
larization observed in SEBN of OMCl not only in silicate is cited in Jena-St.Petersbm-g Database of 

the but also in the L-band, the size distribution Optical Constants (JPDOC, http : //www, astro . uni-jena.de/Labor 

of silicate grains should range from 0.15 to 1.5/um; 
these sizes are larger than those assumed in pre- 
vious papers. With those grains, we deduce the 
Rayleigh reduction factor R > 0.5 in the SEBN 
region. Such a strong alignment cannot be ex- 
plained by the Davis-Greenstein mechanism. We 
suggest alignment by radiative torques as an al- 
ternative mechanism. We also investigate possible 
orientations of grain alignment, and those of in- 
cident and scattered beams in our models. The 
orientations should be almost optimal, and this 
restricts possible configurations significantly. 

3. We investigate how the circular polarization 
in scattered light is affected by linear polarization 
pu in incident light. The effect of linearly polarized 
incident light is small, i.e., the difference between 
circular polarization for linearly polarized incident 
light pii and that for nonpolarized incidence is less 
than Q.ipii in our models. This result shows that 
the conversion from linear to circular polarization 
would not be a dominant process to produce large 
circular polarization, unless pu is extremely large. 
Also the effect of dichroic polarization is small in 
our models, although it will be more significant 
than scattering if the grain sizes are less than those 
assumed here. 

4. If the grains are composed of silicates and 
ice, the effect of the 3/im ice band should appear 
in polarization of scattered light. In our models, 
the degree of linear polarization increases while 
that of circular polarization decreases in the 3/im 
band. Linear and circular polarimetry of the 3/im 
ice band in the SEBN region should be important 
to investigate the details of the scattering process. 
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Fund (FY2008), and the Natural Sciences and 
Research Council of Canada for supporting this 
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A. Sign of the Stokes Parameters 



We use the same definition of polarization as in Ivan de HulstI (|l957l l. i.e., the circular polarization is 
positive {V > 0) when the electric vector is rotating clock wise with time as see n from an observer. However, 
the time component of the electric and magnetic fields in Ivan de Hulsd (|l957( ). i.e., e"'"*'^*, is different from 
that in the FIM and Tmat calculations (e~*"*). The Stokes parameters are thus written for incident light as 



/; — Ei2E*2 + EiiE*i, 

Qi = E.i2E*2 - EiiE*i, 
Ui = Ei2E*i + EiiE*2, 
Vi = -i{EaE*-^ ~ EiiE*^), 



(Al) 

(A2) 
(A3) 
(A4) 



where En and Ei2 ar e amp litudes of orthogonal electric vectors. The sign of Vi in eq.(A4) is different from 
that bv lvan de Hulst (1953) (P-41). The Stokes parameters for scattered light are also written in the same 
manner. The express i ons fo r transformation matrix Fjk of the Stokes parameters (eqU]) are different from 
those in Ivan de Hulst ( 19571 ) (p. 44) for the sign in the fourth line and in the fourth column, except for ^44. 
This definition of the Stokes pararn eter V is different from t hat in our previous papers ([Matsumura fc BastienI 
2004 iBastien fc MatsumuralbOOsh . but is same as that of iGledhill fc McCalll (|2000[ ) and other papers. 



B. Expression of F41 for a Precessing Grain 

When the grain is small, or in the Rayleigh approximation, the component _F4i of the Mueller matrix is 
expressed only with the first term /^^^ of eq.([7]) or eq. pTj) : 

fii2/-'. ^ . X 3/5 



Fii{9i,9s,(j)s) = 0112/ {9i,0s,(j)s) = aii2 • -\/ -sin 6*^ cos 6*^ sin sin (Bl) 

4 V TT 

We now derive the expression of /^^^ when the grain is precessing (S-model) around the direction of _B, and 
prove eg. p^ . We first rewrite eq.(Bl) with the dashed angles of Fig. [31 which are based on the direction of 
alignment B, and then integrate it over (j)'^ from to 27r. The term /^^^ can be written as 

/"^ cxsin6'iCOs6',sinesca(sinZB/S'cosZSM-cosZB/S'sinZSM), (B2) 

where we have used the relation 

sin e, sin 0^ sin Q^ca sin{ABIS - ZBIA), (B3) 
for the spherical triangle IAS. We further use the following relations: 

sin Qsca sin ZBIS = sin 9'^ sin 01, , (B4) 

for the spherical triangle BSI, and 

sin 9i sin ZBIA sin 9'^ sin 4>a , (B5) 

sin 9i cos ZBIA = cos 9'^ sin 9'^ — sin 9'^ cos 6*,- cos , (B6) 
cos 9i = cos 9'^ cos ■ + sin 9'^ sin 6* ■ cos (j)'^, (B7) 

for the spherical triangle BIA. 

With eqs.(B4)-(B7), eq.(B2) can be rewritten as 

/"^oc {A - B cos (Ij'^-C sin (l3'^){D + E cos (j)'^) (B8) 
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where 



A 



cos 9'^ sin 9[ sin 6'^ sin 0^ 



(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 



C 



B 



sin 9'^ cos 9[ sin sin (/)^ , 
sin sin 9sca cos ZBIS, 

D = cos^- cos6'^, 

E = sin6'- sine*;, 



which are constant at present. If we integrate /^^^ over (j)'^ from to 27r, only the constant term and the 
term with cos^ cj)'^ in eq.(B8) remain, and other terms vanish. We thus finally obtain the average /^^^ as 



Comparing eq.(B14) with eq.(Bl), we see the angular dependence is the same for 9'^, 9'^, and (p'^, but the 
amplitude is different by a factor of (1 — (3/2) sin^ 9'^), which is the same as the Rayleigh reduction factor R 
(see eq.®). Therefore, we have shown that < /^^^ > is proportional to R, i.e. ea. p^ in Section r2. 3. 21 

The second term /^^^ in eq.([7]) or eq. (fT7|) is also proportional to the Rayleigh reduction factor. One can 
prove it in the same manner as for /^^^. 



< >cx sin 61^ cos 61^ am9', sin 0^(1 - (3/2) sin^ 9'^). 



(B14) 
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2 4 



Fig. 1. — Comparison of the results computed 
with FIM (symbols) and with Tmat(lines) for an 
oblate particle (axial ratio is 2:1 and m = 1.7). 
The circles (FIM) and the solid line (Tmat) are 
results for (6*^, 6*^, (/"s) = (60°, 60°, 135°), and the 
triangles (FIM) and the broken line (Tmat) are 
those for (6*,, 6*^, 0^) = (60°, 29°, 135°). 



Fig. 2. — Configuration of grain relative to the 
incident and scattered light beams. The grain is 
at the origin and the light source is within the 
plane of y = (and 2 < in this figure). The 
incident light goes in the direction of / within the 
xz-plane and its direction is defined by the angle 0i 
while the scattered light going in the direction of S 
is defined by the angles 9s and (f>s . The scattering 
angle Qsca is the angle between / and S. 




Fig. 3. — Geometry of the problem showing the 
grain configuration with respect to the directions 
of the incident / and scattered S light beams for 
a precessing grain. A is a direction related to the 
symmetry of the grain and B the direction of align- 
ment. See the text for more details. 
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Fig. 4. — Comparison the circula r polarization 
Pc from iGledhill fc McCalll (|2000l ) (circles) and 
that computed here (lines) for spheroidal silicate 
grains (m=1.71+0.03i, axial ratio 2:1). The cal- 
culations were performed at a wavelength of 1 
/Ltm for a size distribution defined by minimum 
and maximum grain sizes of 0.1-1.0 fim and a 
power law with an index of —3.5. For oblate mod- 
els, the incident and scattered directions are set 
to {0'^,e',,(l)',) = (72.0°, 72.9°, 95.7°). For prolate 
models they are {e[,e'„ = (50.0°, 45.6°, 145.3°) 
for open circles and dotted line, and (6*^, 6*1., 0(,) — 
(75.5°, 82.8°, 112.5°) for dashed line. 
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Fig. 5. — Dependence of anmi on \m\xeq or Xeq 
for (a) an oblate grain with an axial ratio 2:1 and 
m = 1.7, (b) same as (a) but for m — 1.3, (c) 
same as (a) but for an ellipsoidal grain with axial 
ratios 2 : : 1, (d) same as (a) but for a prolate 
grain, (e) same as (a) but for a spinning prolate 
grain, (f) same as (a) but for m = 1.1 + O.li, (g) 
same as (a) but for an oblate grain with axial ratio 
3:1, (h) same as (a) but for m = 2 + li, and (i) 
same as (a) but for m = 3 + 3i. Solid line is 
for (n, m, I) — {\, 1, 2), dotted line (2, 1, 1), broken 
Hne (2, 2, 2), long dashed line (2, 1, 3), dotted chain 
(3, 1, 2), and dashed chain (3, 2, 1). 
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Fig. 6. — Dependence of a^T^qi on R. The values 
of (n, m, are (a) (2, 2, 2), (b) (2, 1, 3), (c) (3, 1, 2), 
(d) (3,2,1), (e) (3,2,3), and (f) (3,3,2). In each 
graph, various curves are given for different values 
of n' , m' , I' as indicated. 
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Fig. 7. — The effects of grain precession on cir- 
cular and linear polarization, (a) The means of 
F4i(i?)/F4i(l) for m = 1.7 for 1785 or less direc- 
tions of incidence and scattering (see the text). 
Values of the Raylcigh reduction factor R are 0.25 
(solid), 0.5 (broken) and 0.75 (long dashed) for 
oblate grains (axial ratio 2:1), and -0.5 (dotted 
chain) for prolate grains (axial ratio 2:1) (see the 
text), (b) Same as (a) but for m = 1.3. (c) Same 
as (a) but for the standard deviations, (d) Same 
as (c) but for m = 1.3. (e) The maximum of abso- 
lute values of the circular polarization pc for oblate 
(solid line) and for prolate (broken line) grains, 
and the linear polarization pi at the maximum \pc \ 
for oblate (circles) and prolate (triangles) grains 
with m = 1.7 and R = 1. (f) Same as (e) but for 
m — 1.3. 
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Fig. 8. — The angles for which \pc\ takes its max- 
imum value, (a) The scattering angle Qsca which 
corresponds to maximum values of pc for oblate 
grains with to = 1.7 (circles), oblate grains with 
TO = 1.3 (triangles), spinning prolate grains with 
TO = 1.7 (crosses), and spinning ellipsoidal grains 
with TO = 1.7 (squares), (b) Same as (a) but for 
0i. (c) Same as (a) but for 9s- (d) Same as (a) 
but for (ps- The axial ratios are 2:1 for oblate and 
prolate grains, and 2 : V2 : 1 for ellipsoids. 



Fig. 9. — Linear polarization pi and circular polar- 
ization Pc by oblate vs. the maximum size Umax- 
(a) pi at at A = 2.2/im for axial ratio 2:1 (solid 
line), and 3:1 (dashed line), with R = 1 (circles) 
and R = 0.5 (triangles), (b) Same as (a) but for 
Pc- (c) Same as (a) but for A = 3.6/i-m. (d) Same 
as (c) but for pc- The values of polarization for ax- 
ial ratio 2:1 are averages within {9i, 9s,4>s) = (60 ± 
10°, 60 ± 10°, 135 ± 10°). Those for axial ratio 3:1 
are for (0,, 0,., 0s) = (90± 10°, 50± 10°, 100± 10°). 
These directions are near to those in which pc 
shows its maximum for A = 2.2/Lfm. 
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Fig. 10. — Wavelength Dependence of silicate 
grains with size distribution of aeq = 0.15— 1.5/xm 
for (a) degree of linear polarization pi, (b) posi- 
tion angle PA, (c) degree of circular polarization 
Pc, and (d) ellipticity Pc/pi- The solid and dashed 
lines with open circles are for R = 1.0 and 0.5, re- 
spectively, for oblate grains with an axial ratio 2:1 
and {0i, 9s,(l>s) = (60 ± 10°, 50 ± 10°, 140 ± 10°). 
The broken line with open triangles is for pro- 
late grains with an axial ratio 2:1 for R = 0.5 
and (fl,, 9s,(f>s) = (50 ± 10°, 90 ± 10°, 110 ± 10°). 
The dashed line with open squares is for ellipsoidal 
grains with an axial ratio 2 : \/2 : 1 for i? = 0.5 
and {0„9s, (j)s) = (70 ± 10°, 50 ± 10°, 120 ± 10°). 
The broken line with open diamond is for oblate 
grains with an axial ratio 1.5:1 for R = 0.25 and 
(6*,, 6*^, 0s) = (70 ± 10°, 70 ± 10°, 140 ± 10°). The 
observed values are plotted with filled circles. 



Fig. 11. — The effects of grain shape on polariza- 
tion at A = 2.2/im. (a) The Rayleigh reduction 
factor R necessary to obtain pc = 15%. The open 
symbols are the results of exact calculations, while 
crosses and pluses are those for oblate and prolate 
grains respectively, computed with an approxi- 
mate formula (see the text). The short dashed line 
is derived f rom dichroic extinction a nd th ermal 
emission by iHildebrand &: Dragovan and 
the filled circle is their preferred value R = 0.25. 
(b) The polarization efficiency p/r when the max- 
imum circular polarization pc of 15% is obtained. 
For ellipsoidal grains, the axial ratios are set as 
1 : ^/r : r where r is the ratio of the maximum 
and the minimum radii. The dotted line is the ob- 
served maximum ratio p/r (|jone£>il98i^) . See text 
for more details. 
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Fig. 12. — The effects of polarization {qi, Ui, and 
Vi) in incident light on that of the scattered light 
for silicate grains: oblate with an axial ratio of 
2:1, size distribution of aeq = 0.15 — 1.5/im, the 
Rayleigh reduction factor R = 0.5, and wavelength 
A — 2.2fim. (a) Stokes parameters for linear po- 
larization {q'g and u'^) of scattered light for po- 
larized incidence with {qi,Ui,Vi) = (0.3,0,0) are 
compared with those {qs and Ug) for nonpolarized 
incident light {qi = Ui = Vi = 0). The filled circles 
show the stokes parameter of qs and q'^ , and open 
ones Us and u's- (b) Same as (a) but for circular 
polarization Vs and v's- (c) Same as (a) but for 
((7s, Ms, Vs) = (0, 0.3, 0). (d) Same as (c) but for 
circular polarization Vs and v's- (e) Same as (a) 
but for (gs, Ms, t)s)=(0, 0, 0.3). (f) Same as (e) but 
for circular polarization Vs and v'^. 



Fig. 13. — Contour plots of the number of models 
that satisfy the observed conditions, as a function 
of linear polarization {qi, Ui) in the incident light. 
Silicate grains are assumed: oblate with an axial 
ratio 2:1, size distribution of aeq = 0.15 — 1.5/xm, 
the Rayleigh reduction factor R — 0.7, and wave- 
length A — 2.2/im. (a) Contour plots out of 4913 
models for which pi > 31%. Contour levels are 
800, 1200, 1600, ... and 3200. The darker ar- 
eas show larger numbers, (b) Same as (a) but for 
-5° < PA < 5°. Contour levels are 400, 500, ... 
900. (c) Same as (a) but for pc > 15%. Contour 
levels are 800, 1200, ... 2400. (d) Same as (a) 
but for 31% < pi < 45%, -5° < PA < 5°, and 
12% < Pc < 18%. Contour levels are 5, 15, 55. 
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Fig. 14. — The ice band feature for oblate 
grains with an axial ratio of 2:1 for (a) degree 
of linear polarization pi, (b) position angle PA, 
(c) degree of circular polarization pc, and (d) flux 
normalized at 2/im. The solid line is for grains 
composed of bare silicate with a size distribu- 
tion of 0.15 — 1.5^m, R = 0.5, and {9i,9s,(f>s) = 
(60 ± 10°, 50 ± 10°, 130 ±10°). The short dashed 
line is the same but for a silicate-ice mixture whose 
volume of ice is 10%. The long dashed line is the 
same but for an ice volume of 25%. The dotted 
chain is the same but for an ice volume of 50%. 
The double dotted chain is for pure ice grains with 
size distribution of 0.193 — 1.93/im, R ~ 0.7, and 
{e^,es, cj)s) = (75 ± 10°, 70 ± 10°, 140 ± 10°). The 
observed values are plotted with filled circles. 
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Table 1 
Conditions for or Pc = 



0° or 180° 
any 
any 
90° 



any 
and 

0° or 180° 
90° or any" 



any 
0° or 180° 

any 
any 



'^For the Rayleigh scattering, 6a can be 
any. 
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Table 2 

Models that Explain Observations in the K-band" 



R 




N 


9' 


"a 




PA' 




Oblate 2:1, L 














0.2 


0.07 


Q 












0.4 


0.15 


14 


60( 8) 


56( 6) 


136( 8) 


46( 7) 


104( 9) 


0.6 


0.24 


23 


63(12) 


48(11) 


140(12) 


49(10) 


102(15) 


0.8 


0.33 


21 


67(15) 


41(14) 


144(13) 


50(12) 


101(19) 


1.0 


0.43 


11 


70(29) 


51(27) 


152(12) 


55(16) 


115(36) 


Oblate 3:1, L 














0.4 


0.12 





- 


- 


- 


- 


- 


0.6 


0.19 


31 


80(16) 


52 ( 7) 


128(11) 


27(13) 


112(14) 


0.8 


0.28 


9 


83(22) 


43( 9) 


135(13) 


35(15) 


113(21) 


1.0 


0.41 


7 


IOC/ 
125(20) 


32( 8) 


114(27) 


-7(26) 


132(21) 


Oblate 2:1, S 














0.4 


0.12 





- 


- 




- 




0.6 


0.19 


31 


72(12) 


49 ( 7) 


1 O/i /I ^^ 
124\^iU ) 


31(10) 


1 no/ 1 o^ 
iU2^i2 J 


0.8 


0.26 


12 


68(17) 


41 ( 7) 


1 Q^!/ 1 o^ 
loD(^i2 J 


43(11) 


99(17) 


1.0 


0.33 


10 


67(15) 


38 ( 6) 


1 "^fi?! 


46(12) 


97(15) 


Prolate 


2:1, L 














0.4 


0.11 















0.6 


0.16 


8 


63(10) 


56( 9) 


128( 6) 


38( 7) 


102(11) 


0.8 


0.23 


14 


59(16) 


61(16) 


121( 8) 


36(11) 


98(16) 


1.0 


0.35 


20 


67(26) 


60(26) 


110(11) 


24(19) 


94(21) 


Prolate 


3:1, L 














0.6 


0.08 















0.8 


0.14 


10 


72( 8) 


72( 8) 


90( 0) 


12( 8) 


84( 3) 


1.0 


0.33 


18 


78(12) 


75(15) 


94(13) 


7(11) 


91(13) 


Prolate 


2:1, S 














0.6 


0.11 















0.8 


0.15 


47 


69(11) 


66(10) 


107( 9) 


20(10) 


96(10) 


1.0 


0.21 


83 


75(16) 


62(16) 


108(12) 


16(14) 


98(14) 


Ellipsoid 2 : %/2 : 


1, L 












0.4 


0.14 















0.6 


0.22 


5 


49( 8) 


56(10) 


140( 6) 


54( 5) 


96(11) 


0.8 


0.31 


8 


69(14) 


38(10) 


132(12) 


40(11) 


96(15) 


1.0 


0.41 


21 


59(27) 


57(32) 


131(13) 


42(18) 


101(32) 


Ellipsoid 2 : \/2 : 


1, S 












0.4 


0.10 















0.6 


0.16 


40 


70(11) 


53( 8) 


120(10) 


28( 9) 


100(11) 


0.8 


0.23 


29 


70(14) 


42( 9) 


127(12) 


35(11) 


97(14) 


1.0 


0.30 


21 


75(10) 


37( 6) 


128(12) 


33(11) 


98(12) 



^The letter 'L' stands for the size distribution of agq = 0.15 — l.Sjum, and 
the 'S' for Ogq = 0.1 — ljum, respectively. The observed values in the K-band 
(A = 2.2iJ,Ta) for p^''^, p^^^ , and PA°^^ are set as 15%, 38%, and 0°, and the 
ranges of acceptance specified by 5p'^^ , ^Pi^^ and 5PA°^^ as 3%, 7%, and 
5°, respectively. The ranges are larger than the observational errors. 
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Table 3 

Models that Explain Observations in the L-band" 



R 




N 


e'. 






PA' 




Oblate 


2:1, L 














0.2 


0.05 





- 


- 


- 


- 


- 


0.4 


0.10 


71 


78(13) 


48(11) 


117(11) 


21(11) 


101(13) 


0.6 


0.14 


62 


77(18) 


40(11) 


129(17) 


33(16) 


102(19) 


0.8 


0.19 


48 


76(24) 


39(12) 


137(20) 


41('20') 


105(25) 


1.0 


0.25 


44 


76(35) 


42(17) 


135(31) 


39(30) 


106(36) 


Oblate 


3:1, L 














0.2 


0.06 





_ 


_ 


_ 


_ 


_ 


0.4 


0.11 


73 


87(17) 


51(11) 


115(11) 


13(13) 


107(14) 


0.6 


0.16 


66 


84(22) 


44(11) 


123(15) 


23(16) 


108(21) 


0.8 


0.20 


40 


94(30) 


38( 9) 


122(27) 


17(27) 


113(30) 


1.0 


0.24 


16 


100(29) 


37( 6) 


119(39) 


12(36) 


115(32) 


Oblate 


2:1, S 














0.4 


0.06 





- 


- 


- 


- 


- 


0.6 


0.09 


7 


112( 9) 


61( 9) 


111( 4) 


-16(10) 


118( 5) 


0.8 


0.12 


12 


119(15) 


47( 9) 


106( 9) 


-18(13) 


120(11) 


1.0 


0.15 


15 


125(15) 


39(10) 


98(17) 


-24(15) 


121(15) 


Prolate 2:1, L 














0.4 


0.05 















0.6 


0.07 


20 


68( 9) 


57( 9) 


110( 7) 


24( 7) 


93( 9) 


0.8 


0.10 


75 


73(16) 


59(13) 


109(11) 


18(13) 


97(14) 


1.0 


0.15 


81 


83(19) 


48(13) 


103(17) 


8(15) 


95(18) 


Prolate 3:1, L 














0.6 


0.05 















0.8 


0.09 


42 


80(10) 


60(11) 


82( 8) 


0(9) 


78( 9) 


1.0 


0.16 


56 


89(12) 


56(14) 


80(16) 


-10(13) 


82(15) 


Prolate 2:1, S 














1.0 


0.07 















Ellipsoid 2 : V2 : 


1, L 












0.2 


0.04 















0.4 


0.08 


13 


78(12) 


46( 6) 


118(15) 


23(13) 


100(13) 


0.6 


0.12 


39 


79(16) 


38(11) 


118(21) 


23(18) 


98(18) 


0.8 


0.17 


34 


75(23) 


38(15) 


127(21) 


33(19) 


99(24) 


1.0 


0.21 


32 


78(31) 


40(18) 


117(37) 


22(32) 


97(34) 


Ellipsoid 2 : V2 : 


1, S 












0.6 


0.07 















0.8 


0.09 


1 


129( 0) 


51 ( 0) 


101( 0) 


-32( 0) 


121( 0) 


1.0 


0.11 


4 


135(13) 


41( 0) 


93(11) 


-37(10) 


123(10) 



'^Thc letter 'L' stands for the size distribution of a^g — 0.15 — and 
the "S" for a^q — 0.1 — respectively. The observed values in the L-band 

(A ^ 3.6/im) for p'^^^ , pf^^, and PA""^^ are set as 9.4%, 57%, and 0°, and the 
ranges of acceptance specified by 5p°^^ , Sp^^^ , and 5PA°^^ as 3%, 11%, and 
5°, respectively. Ths raqnges are larger than the observational errors. 
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